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THE  BULIKEN  CALIBRATia;  SHOT  IN  SOUTHEAST  MISSOURI 


Abstract 


On  June  28,  1963,  a  20,000  lb.  high-explosive  ’’calibra¬ 
tion-shot,”  code  named  BILLIKEN,  was  detonated  near  the 
epicenter  of  the  magnitude  5  Southeast  Missouri  earthquake  of 
March  3,  1963.  The  experiment  was  suggested  by  personnel  of 
Stanford  Research  Institute  after  field  units  set  up  in  the 
epicentral  area  failed  to  record  aftershocks.  The  shot  was 
sponsored  by  the  Advanced  Research  Projects  Agency  and  con¬ 
ducted  jointly  by  the  Crustal  Stxidies  Branch  of  the  U,  S. 
Geological  Survey  and  Saint  Louis  University. 

Using  arrival  times  at  stations  co!rparable  to  those  used 
in  locating  the  March  3  earthquake,  the  instrumental  epicenter 
of  BILLIKEN  was  located  within  2  to  U  kilometers  of  the  shot 
point.  This  confirms  the  accuracy  of  epicenter  locations  in 
this  area  from  the  data  of  near-regional  stations.  The  shot 
has  also  permitted  the  "calibration”  of  the  area,  making  pos¬ 
sible  epicenter  locations  using  time  differences  of  arrivals 
at  pairs  of  stations.  To  the  east,  north  and  west  of  the 
shot-point  travel  times  from  BILLIKEN  agree  with  those  of  pre¬ 
vious  earthquake  studies  and  of  a  current  refractions  survey; 
to  the  south  negative  residuals  were  obtained.  A  comparison 
of  the  (Pg-Pn)  vs  (Sg-Pg)  tii.e  intervals  from  BILLIKEN  with 
those  for  the  earthquakes  of  2  February  1962  and  3  March  1963 
indicates  a  focal  depth  deep  in  the  crust  for  the  former, 
shallow  for  the  latter. 


Ihe  BILL13CM  Calibration  Shot  in  Southeast  Missouri 

by  Williiac  itauderf  S»  John  Dowlin^^ 
and  Wayne  Jaclason 


!•  Introduction 

March  3*  1963,  a  magnitude  3  earthquake  occurred  in  the  southeast 
Missouri  area.  Within  hours  of  the  occurrence  of  the  earthquake  word  was 
sent  to  Stanford  Research  Institute  giving  the  approximate  magnitude  and 
location  of  the  shock.  Immediately,  as  a  part  of  its  study  of  after¬ 
shocks  of  earthquakes,  Stanford  Research  Institute  took  steps  to  estab¬ 
lish  tMporary  monitoring  stations  in  the  epicentral  area*  Three  such 
stations  were  established  and  monitoring  of  aftershocks  began  on  the 
fifth  day  after  the  ir.'ia  shock.  Within  the  first  twenty  days,  however, 
no  aftershocks  were  recorded  (see  Lange  and  Westphal,  1963)*  The  failure 
to  record  aftershocks  occasioned  the  eaqpeJlment  which  is  here  reported. 

It  was  reasoned  that  the  failure  to  record  aftershocks  could  be 
explained  by  one  of  two  cr  ^es:  1)  no  aftershocks  actually  occurred; 

2)  the  monitoring  teams,  due  to  errors  in  the  location  of  the  epicenter 
of  the  earthquake,  were  not  looking  in  the  correct  place.  Accordingly, 
in  a  letter  of  March  20,  1963,  to  the  Advanced  Research  Injects  Agency, 
Stanford  Research  Institute  suggested  that  the  circumstances  in  south¬ 
east  Missouri  afforded  good  opportunity  for  a  chemical  hi^-exploeive 
oU.ibrati^  shot.  The  purpose  of  the  experiment  would  be  to  confirm 
the  accuracy  of  the  location  of  the  epicenter  of  the  earthquake  of  March 
3,  1963,  and  to  make  possible  a  refinement  of  the  location  of  the 
epicenters  of  future  earthquakes  which  mi^t  occur  in  the  same  region. 

It  was  also  suggested  that  the  calibration  shot  would  provide  an  ero¬ 
sion  "seismic  signature"  to  serve  as  a  standard  for  the  region  or  as  a 
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criterion  to  distinguish  a  given  event  as  an  earthquake  or  as  an  explo¬ 
sion. 

Following  the  initial  suggestion  by  Stanford  Research  Institute «  on 
April  9  the  Advanced  Research  Projects  Agency  sent  notice  of  the  plan  to 
^e  Crustal  Studies  Branch  of  the  U.  S.  Geological  Survey,  which  was 
scheduled  to  conduct  a  cooperative  crustal  refraction  survey  in  Missouri 
with  Saint  Louis  University  the  latter  part  of  June.  ARPA  suggested  to 
the  USQS  that  it  ni^t  conduct  the  calibration  experiment  at  tnls  same 
time  and  requested  a  proposal  to  this  Intent.  A  proposal  was  submitted 
on  May  l6th,  the  experiment  was  authorized  a  few  days  later,  and  the  cal¬ 
ibration  shot  was  scheduled  for  June  28th.  13ie  exiJ.08ion  took  place  on 
schedule  at  5  A.M.  £ST,  on  that  date.  The  shot  was  code  named  BILLIKIlf. 

2.  Shot  Point  and  Recording  Operations 
The  epicenter  and  the  isoseismal  lines  of  the  earthquake  of  March 
3 1  1963  are  indicated  in  figure  1.  Ihe  earthquake  was  felt  quite  widely, 
the  radius  of  the  felt  region  being  about  200  miles.  As  indicated,  the 
intensity  in  the  immediate  epicentral  region  was  VII  on  the  Modified 
Mezcalli  scale.  The  region  of  the  epicenter  is  the  alluviated  bottom¬ 
lands  lying  midway  between  Poplar  Bluff  and  New  Madrid,  Missouri.  The 
calibration  shot  was  planned  to  be  located  as  nearly  as  possible  at  the 
epicenter  of  the  earthquake. 


Recording  Plan 

Several  government  organizations,  universities,  and  private  corpora¬ 
tions  collaborated  in  the  recording  of  the  event.  The  plan  for  the  ve- 
cording  is  illustrated  in  figure  2.  Since  the  explosion  was  intended  as 
a  check  on  the  accxiracy  of  the  epicenter  location,  this  aspect  of  the 
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experiment  could  best  be  obtained  by  recording  at  the  permanent  seismo- 
graphic  stations  in  the  region.  These  are  indicated  by  the  black  dots 
on  the  map.  With  the  exception  of  Cape  Girardeau,  the  nearest  seismo- 
graphic  station  is  at  an  epicentral  distance  of  over  200  kilcmeters.  The 
magnification  at  most  of  these  stations  is  of  the  order  of  50«000.  Con¬ 
sequently,  it  was  feared  that  the  signals  of  the  explosion  mi^t  not  have 
sufficient  amplitude  to  be  recorded  at  the  permanent  stations.  This  cir¬ 
cumstance  controlled  one  aspect  of  the  dispersal  of  the  available  recording 
units;  certain  of  the  more  sensitive  field  recording  units  were  located 
either  side  of  or  near  the  permanent  stations.  It  also  sea.-ed  advisable, 
as  much  as  possible,  to  establish  continuous  lines  in  x:rious  directions 
from  the  shot  point;  tliis  consideration  provided  a  second  control  over 
the  disposal  of  the  recording  units. 

The  foremost  and  hipest  quality  field-recording  potential  was  that 
provided  by  the  preser-e  of  the  field  teams  of  the  U.  S.  Geological  Sur¬ 
vey  Crustal  Studies  Branch.  Ten  recording  trucks  were  in  the  region  for 
the  crustal  refraction  survey  in  Missouri.  Each  truck  was  capable  of 
putting  out  six  vertic«tl  seismometers  at  500  meter  intervals  along  a  line 
2.5  kilometers  in  length.  These  USGS  units  were  selected  for  dispersal 
so  as  to  provide  travel  time  information  in  the  vicinity  of  the  permanent 
seismographic  stations  of  the  Saint  Louis  University  net.  At  Holla  and 
Bloomington  one  USGS  field  unit  was  placed  about  20  kilometers  either 
side  of  the  permanent  station.  At  stations  more  distant  from  the  shot 
point  the  field  units  were  located  along  lines  in  the  direction  of  the 
station  fran  the  shot  point. 

In  addition  to  the  penaanvit  stations,  Saint  Louis  University 
established  four  temporary  recording  units  with  Beniof*  seismometers  and 


Figure  2,  Perinanent  stations  and  disposition  of 
recording  tears  for  the  EILLIKiJi  calibration 
shot  of  28  June  1963. 
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with  conventional  photographic-recording  used  in  earthquake  seismology. 
Three  of  these  units  were  located  along  a  line  between  St.  Louis  and  the 
shot  point.  The  fourth  was  located  in  sourthem  Illinois  along  a  line 
towards  Bloomington,  Indiana. 

Two  field  units  were  made  available  by  the  U.  S.  Coast  and  Geodetic 
Survey.  These  were  located  at  60  and  85  km  from  the  shot  point  in  the 
direction  towards  Holla.  At  these  sites  the  recording  spread  consisted 
of  six  two  cycle  per  second  seismometers  located  at  the  comers  of  a 
hexagon  whose  lateral  extent  was  300  feet. 

The  Geotechnical  Corporation  fielded  two  single  component  units. 

The  first  of  these  occupied  a  former  LRSM  site  near  Parsons,  Tennessee; 
the  second  was  located  near  Huntsville,  Arkansas.  The  University  of 
Michigan  provided  one  field  station  vdaich  was  located  near  Corinth, 
Mississippi.  The  Geological  Surveys  of  Illinois  and  Indiana  each 
fielded  a  seismic  refj*action  unit;  these  were  located  near  C3iampaign- 
Urbana  and  near  Mt.  CauTnel,  Indiana,  respectively. 

Thus,  as  much  as  possible,  there  was  provided  good  recording  at 
or  near  the  permanent  seismographic  stations  of  the  region,  and  recording 
along  lines  extending  in  several  directions  from  the  shot  point. 

Shot  Point 

The  RTT.T.TKCT  shot  point  was  located  about  three  km  from  the 
ins trumen tally  located  epicenter  of  the  March  3rd  earthquake.  The 
explosion  consisted  of  a  total  of  l8,400  pounds  of  DuPont  Super 
’’Tovex"  gel  fired  in  a  10-i»«ik  hole  pattern.  The  shot-hole  pattern  (see 
figure  3)  was  designed  to  provide  some  cancellation  of  the  horizontal 
components  of  surface  waves  to  minimize  destructive  effects  of  the  ex¬ 
plosion  on  nearby  structures.  Six- inch  holes  were  drilled  to  a  d^pth  of 


Holes  A  thru  I  •  2000  lbs  -r  primers 
Hole  J  •  400lbs  -t  primer 

Depth*  115-120'  center  of  sperical  charge 


Plan  of  the  shot-hole  pattern  at  ti 
BILLIKEIJ  shotpoint. 
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120  ft  and  their  bottoms  were  enlarged  by  prolonged  washing  with  the 
drill  bit.  The  material  drilled  through  was  primarily  saturated  clays  on 
the  surface  trending  toward  saidy  clay  to  pure  sand  at  the  base  of  the 
holes.  Super  "Tovex,"  which  is  a  hi^-energy,  high-density  explosive, 
was  bulk  loaded  in  order  to  fill  the  entire  cross-section  of  the  hole. 
Holes  A  throu^  I  (figure  3)  were  loaded  with  2,000  pounds  each  and  hole 

J  was  loaded  with  400  pounds.  The  density  of  Super  "Tovex"  is  1.7 

3  3 

gm/cm""^  and  the  wei^t  is  approximately  100  pounds/ft"^. 

The  charge  was  detonated  automatically  by  a  firing  circuit  activated 

by  an  electronic  chronometer  phased  to  coincide  closely  with  WWV  timing 

signals.  The  time  break,  radio  tones  and  other  instrumental  and  timing 

information  were  recorded  on  photograjiiic  paper  operating  at  5  inches/sec. 

The  time  of  detonation  can  be  determined  to  within  5  millisec  or  less. 

The  shot  point  coordinates  including  location  and  shot  time  are  as 

follows:  Latitude  36®  43.00*  N 

Longitude  90®  07.62'  W 
Shot  time  10-00-00.13  GOT 

3.  Interpretation  of  Results 

Signals  were  recorded  from  the  BILLIKEN  calibration  shot  at  29  sta¬ 
tions,  The  time  of  arrival  of  the  first  recorded  P  motion  is  given  in 
Table  1,  The  times  reported  at  stations  occupied  by  the  USGS  field  units 
are  the  arrival  times  at  the  first  and  the  last  seismometer  of  the  spread. 
At  several  stations,  e.g,,  Mt.  Sterling,  Ulinois,  did  not  record;  the 
first  recorded  motion  at  these  stations  is  listed  in  the  column  headed 

>  ",  The  residuals  ("M")  given  in  the  last  column  of  the  Table 
S 

correspond  to  the  difference  between  the  observed  arrival  time  and  the 
arrival  time  computed  by  the  travel  time  oirves  used  in  the  epicenter 


Table  1,  Recording  sites  and  time  of  arrival  of  first  recorded 
P  motion,  BILLDOJ  shot,  28  June  196 >. 


Station 

Group 

P  ,  P*,  P 

g’  n 

P 

At 

misinor.  Mo. 

USCGS 

59.32 

10-00-10.16 

+0.25 

Cape  Girardeau,  Ho. 

SLU 

85.00 

10-00-14.10 

■fO.03 

Ellington,  Mo. 

USCGS 

88.02 

10-00-1^.41 

-0.13 

Fredericktovm,  Mo. 

SLU 

94.85 

10-00-16.10 

+0.45 

Flat  River,  Mo. 

SLU 

129.57 

io-oo-a.8o 

*0.74 

DeSoto,  Mo. 

SLU 

161.39 

10-00-25.50 

-0.75 

Salem,  Mo. 

uses 

168.06 

10-00-27.62 

+0.30 

Salem,  Mo. 

USGS 

170.42 

10-00-27.95 

Berry  Farm,  HI. 

SLU 

i88.a 

10-00-30.70 

+0.25 

Tyson,  Mo. 

SPRENG 

203.41 

10-00-32.70 

+0.10 

Rolls,  Mo. 

SLU 

203.98 

10-00-33.10 

+0.50 

Parsons,  Tenn. 

GEO 

209.46 

10-00-31.40 

-2.30 

St.  Louis,  Mo. 

SLU 

a3.26 

10-00-35.00 

+1.10 

Searcy,  Ark. 

USGS 

221.88 

10-00-36.47 

+1.57 

Vienna,  Mo. 

USGS 

220.77 

10-00-35.10 

+  .a 

Vienna,  Mo. 

USGS 

223.25 

10-00-35.35 

Florissant,  Mo. 

SLU 

232.40 

10-00-37.50 

+0.30 

Corinth,  Miss. 

MICH 

267.18 

10-00-38.20 

-2.30 

Mt.  Carmel,  Ind. 

nnxss 

287.45 

10-00-45.07 

+2.0 

Little  Rock,  Ark. 

SLU 

294.00 

Huntsville,  \zic. 

GEO 

326.68 

10-00-49.80 

+2.0 

Mt.  Sterling,  HI. 

US(^ 

350.27 

10-00-54.01 

+3.41 

Mt.  Sterling,  HI. 

USGS 

352.72 

10-00-54.30 

Loogootee,  Ind. 

USGS 

354.85 

10-00-52.36 

+1.06 

Loogootee,  Ind. 

USGS 

356.58 

10-00-52.58 

Fayetteville,  Ark. 

UARK 

370.58 

10-00-01.2 

+8.2 

(Hinton,  Mo. 

USGS 

373.56 

10-00-54.93 

+1.13 

(Hinton,  Mo. 

USGS 

375.76 

10-00-55.18 

Bloomington,  Ind. 

SLU 

420.31 

10-01-06.27 

+7.10 

McMinnville,  lYi. 

CPSO 

428.43 

10-(H-01.00 

0 

. 

+ 

Roseville ,  Ill . 

USGS 

443.76 

10-01-02.97 

+  .77 

Roseville,  HI. 

USGS 

446.  a 

10-01-03.16 

Franklin,  Ind. 

USGS 

463.38 

10-01-05.63 

+1.13 

Franklin,  Ind. 

USGS 

464.98 

10-<H-05.75 

Olathe,  Kans. 

USGS 

471.15 

10-01-05.58 

+0.08 

(Ha the,  Kans. 

USCiS 

472.56 

10-01-05.81 

Wartburg,  Tn. 

LRSM 

486.a 

10-01-17.50 

+10.2 

Manhattan,  Kans. 

SLU 

630.38 

Dubuque,  la. 

SLU 

644.62 

Grapevine,  Tx. 

GEO 

758.62 

10-01-41.30 

+0.7 

DESOTO  28  JUNE  1963 
D  =  161  KM 


Z  STRIP  MINE  EXPLOSION,  D  =  68km 


Z  LOCAL  EARTHQUAKE,  D  =  75  km  f 


Figure  It.  Copparison  of  records  from  BILLIKciN  and 

ot/her  small  seisniic  sources  of  the  Missouri  area 
Ua  (top)  BHiLIKEti  record  at  deScto,  Missouri, 
hb  (mdddle)  strip  nine  explosion,  also  recorded 
at  DeSoto,  Uc  (bottom)  small  local  earthquake 
recorded  near  Poplar  B3uff, 
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determination.  'Hie  records  obtained  were  sufficient  to  achieve  the 
principal  objectives  of  the  experiment.  Particularly  at  the  larger 
distances,  however,  the  recordings  are  minimal  and  the  data  are  of 
limited  usefulness  for  certain  types  of  studies. 

Lliaracter  of  the  Seismic  Signature 

Figure  4a  shows  the  recording  of  the  BFilKEN  event  at  one  of  the 

temporary  stations,  using  traditional  earthquake  recording  equipment. 

The  record  shown  was  written  at  DeSoto  at  a  distance  of  l6l  kra  from 

the  source.  The  paper  speed  was  1?0  in  per  minute,  dear  P  and  S 

S  S 

phases  are  here  recorded.  It  is  noteworthy  that  no  significant  sur¬ 
face  wave  train  appears  on  the  record.  This  is  contrary  to  the  signa¬ 
tures  nomally  obtained  from  explosions  in  this  region.  Figure  4b,  for 
instance,  shows  the  characteristic  record  which  is  obtained  from  strip 
mine  explosions  in  the  area,  whereas  figure  4c  shows  a  record  which  is 
characteristic  of  small  local  eeirthquakes .  As  a  general  rule  records 
from  industrial  explosions  in  the  Illinois -Missouri  aurea  feature  a  well- 
developed  short  period  surface  wave  train,  sucn  as  that  seen  in  figure 
4b,  which  may  eeisil”-  '^e  recorded  to  distances  of  200  km  and  greater. 

The  records  from  small  earthquakes,  on  the  other  hand,  are  notably  lacking 
in  the  surface  wave  trains.  The  contrast  between  the  BTTJiIKSN  explosion 
record  and  those  from  industrial  explosions  may  well  be  due  to  the  source 
environment.  BILLIKEN  was  detonated  in  alluvium  overlying  recent  loosely 
consolidated  sediments.  The  industrial  explosions  normally  occur  in  well 
indurated  carbonate  rock. 

Epicenter  Location 

Signals  identified  as  P  wave  first  arrivals  were  recorded  from 
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BILLIKI2I  at  22  stations.  In  order  to  achieve  the  principal  purpose  of 
the  experiment  the£?»^  data  were  used  in  an  epicenter  location  program. 
Figure  5  will  illustrate  this  phase  of  the  data  interpretation.  The 
cross  on  the  figure  marks  the  location  of  the  instrumental  epicenter 
for  the  earthquake  of  March  5i  19^5 i  as  located  independently  by  the 
U,  S,  Coast  and  Geodetic  Survey  and  by  Saint  Louis  University.  Each 
of  these  determinations  was  to  the  nearest  tenth  of  a  degree  and  the 
two  determinations  were  in  exact  agreement.  The  BILLIKEN  shot  location, 
as  indicated  on  the  figure,  was  situated  as  close  to  the  instrumental 
epicenter  as  was  practical, 

A  second  instrumental  epicenter  of  the  March  3rd  earthquake  was 
determined  by  Saint  Louis  University  after  a  computer  epicenter  determin¬ 
ation  program  was  developed  late  in  the  summer  of  1963»  Using  data  from 
seven  stations  -  five  stations  of  the  Saint  Louis  University  net  (St. 
Louis,  Florissant,  Manhattan,  Little  Rock,  HLoaaington)  plus  the  arrival 
time  at  Lawrence  and  at  the  Cumberland  Plateau  Seismological  Observatory 
(CPSO),  McMinnville,  Tennessee  -  a  revised  epicenter  was  determined.  In 

this  determination  the  travel  time  curves 

P  =  +  ^  ,  0-iA<150 

P*=3.2  +  L,  ,  130<.C^<170 

P^=6.0  +  V8.24  ,  6  >170 

were  used.  These  travel  time  curves  correspond  to  a  focal  depth  of  15- 

20  km.  The  least  squares  program  as  used  in  this  determination  holds 

the  focal  depth  fixed.  The  resulting  epicenter  was 

36.649‘-N,  90.026®W, 

0  =  17-50-12.66  GCT 

Std,  dev,  of  travel  times  =  ,l4  sec. 

This  epicenter  is  indicated  on  the  map  in  figure  5  Uy  the  black  dot.  This 
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90.20*  90.15*  90.10*  90.05*  90.00* 


36.75* 


36.70* 


36.65* 


36.60* 


Figure  5«  Index  map  of  southeast  Kissouri  area  showing 
the  location  of  epicenters  of  the  earthquake  of 
3  Karch  1963  and  of  BILLIKcI^’  (see  text). 


Table  2.  Tabulation  of  data  relevant  to  the  determination 


01 


the  epicenter  of  the 


BILLIK£,N 


event. 


ST  LOUIS  UNIV  EPICENTER 

PROGRAM,  AUG 

1963-  TVM,  MLG 

DATE- 

06286S  ASSUMED 

DEPTH  0 

- 

“vi  f\WI^5€C 

11  km 

TRAVEL  TIME  DATA  USED  AT 

LESS  THAN  1000 

KM 

VEL,  1 

<M/SEC  6.10  6.40 

6.70  8.15 

TIME 

intercepts  0.00 

.92  2.27 

7,57 

XtrlVTn 

-“OrrtvnvScC 

CROSS 

OVER  DIST,  KM  l?0. 

00  193.00  200 

1* 

- 

b.iD  twvjKr 

EPICENTER  TIME 

LATITUDE 

STD,  DFV. 

NO.  TRIALS 

BILL 

10  00  00.10 

36.6821 

90.108 

,87562 

ST  A 

DIST. -KM  D 

1 ST,-DFG 

ARR  TIME 

TRV  TIME-SFC 

DFV-SFC 

ELLS 

62.2936 

.560? 

1 

10.16 

10.05 

-.1608 

COM 

87.078? 

.7831 

1 

14.10 

13,98 

-.2855 

ELLI 

91.7322 

.8250 

1 

14,41 

14.29 

-.7388 

rPED 

98.9i9ft  " 

.8W7 

1 

~  16',  1ft 

l5,99 

“ *2266 

FLAT 

133.7705 

1.2030 

1 

21,80 

21.69 

-.1316 

sTlT 

T 

27.62 

2  7,51 

—.2565 

BERR 

189.2798 

1.7023 

1 

30.70 

30.58 

.0942 

TYSN 

207.5398 

1,8665 

1 

32.70 

32.59 

-.4442 

ROL 

207.8233 

1.8690 

1 

33.10 

32.99 

-.0794 

SLM 

217.1895 

1.9533 

1 

35.00 

34.89 

.6721 

SEAR 

220.1802 

1.980? 

1 

36.47 

36.36 

1,7746 

VIEN 

227.1634 

2.0430 

1 

35.10 

34,99 

-.4508 

FLO 

236.3488 

2.1256 

1 

37,50 

37,38 

.8197 

LOOG 

357.7710 

3.2176 

1 

52.36 

52,25 

.7821 

CLIN 

377.0583 

3.3911 

1 

54,93 

54,82 

.9859 

CPSO 

42T.T072 

3.'B2T? 

T~ 

“i  01  .oO 

"66,88 

- 1  .0991 - 

ROSE 

447.8199 

4.0275 

1 

1  02.97 

62,85 

.3423 

FRAN 

464  66iT5 

4.170^ 

1 

65.51 

- rrni - 

OLAT 

474.6101 

4,2684 

1 

1  05,58 

65.46 

-.334,7 

OESO 

165.5544 

1,4889 

1 

25,50 

25,39 

-1.3970 

PARS 

206.0667 

1 ,8532 

1 

31,40 

31.29 

-1.563C 

CORI 

263,2134 

2.3672 

1 

38.20 

38,09 

-1.7745 

revised  epicenter  is  approximately  the  center  of  gravity  of  several 
determinations  at  oouthem  Methodist  University,  (Tag^u't,  personal  communi¬ 
cation)  •  It  is  distant  about  8  km  fiw  *he  first  epicenter  location. 

In  order  to  locate  the  BILLIKi}!  event  after  the  manner  in  which  an 
earthquake  ^n  this  regicm  would  be  located,  the  arrival  times  were  selected 
from  five  stations  (RoUa,  St.  Loviis,  Florissant,  Searcy,  Franklin;  also 
CFSO)  of  the  Saint  Louis  University  network,  or,  where  one  or  other  of 
these  stations  failed  to  record  tiie  event,  from  one  of  the  nearby  USGS  field 
units,  ihe  arrival  time  from  McMinnville  was  also  used.  Fran  these  data, 
using  the  same  travel  time  curVw>8  as  in  the  location  ox  the  revised 

earthquake  epicenter,  the  following  location  for  RTT.T.Tktm  was  obtained: 

36.725®  N,  90.108®  W 

0  =  05-00-03.32 

Std.  dev.  of  travel  time  s  0.l4  sec. 

This  location  is  indicated  oy  the  inverted  trian^e  on  the  map.  Using  a 
sli^tly  different  velocity  ox'  8.10  km/sec.  and  a  time  intercept  to 
correspond  to  a  surface  focus,  data  from  these  same  stations  gave  the 
second  Saint  Louis  University  epicenter  indi  ;ed  on  the  figure.  These 
epicenters  are  2  and  4  km  distact  from  the  shot  point,  respectively. 

To  make  better  xusf.  of  all  of  the  data  from  the  BILLIKEN  experiment 
a  third  epicenter  determination  was  made  using  the  data  fron  all  22  sta¬ 
tions  which  recorded  P  wave  first  arrivals.  These  data  are  tabulated  in 
Table  2.  The  previous  epicenters  were  located  using  a  two-layer  crust 
with  a  6,08  and  6.4l  velocities,  resp'ccively.  As  a  result  of  the 
crustal  refraction  survey  cwiducted  during  the  suimner  of  1963,  and  also 
as  a  result  from  surface  wave  studies  in  the  region,  a  better  knowledge 
of  a  crustal  structure  is  now  available.  The  more  recent  studies 
indicate  the  probable  presence  of  velocity  gradients  within  the  crust. 
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or,  if  the  crxistal  structure  is  to  be  approximated  by  constant  velocity 
layers,  that  a  three-layer  crust  is  probably  more  applicable*  Accordingly, 
a  three-layered  crust  with  velocities  as  indicated  in  Table  2  was  used  in 
determining  the  epicenter  frcwj  the  P  arrival  times  listed  in  the  table. 

The  resulting  epicenter,  given  by  the  open  circle  nximbered  1  in  the 
figure  is 

36.682®  N,  90.108®  W 
0  =  10-00-00.10 

Std,  dev.  of  travel  times  =  O.876  sec. 

The  last  three  stations  listed  in  the  table  have  larg'-r  than  average 
residuals  in  their  travel  time  from  the  actual  shot  point  and  shot  time 
(see  below)  than  do  the  other  stations.  Rejecting  these  last  three 
stations,  another  epicenter  determination  was  made  using  the  same 
crustal  model.  This  gave  the  epicenter,  marked  by  the  open  circle 
niunbered  2 

36.7^3*  N,  90.135*  W 
0  =  10-00-01.00 
Std.  dev,  =  0.460  sec. 

All  of  tfie  above  epicenters  fall  v/ithin  2  to  4  km  of  the  actual 
BILLIKEN  shot  point.  As  a  result  it  is  concluded  th.\t  the  reliability 
of  the  travel  time  curves  which  have  been  developed  for  this  region  is 
established,  and  that  by  using  these  curves  epicenter  determinations  can 
be  made  from  au:rival  times  at  local  and  regionaJ.  stations  with  an  accuracy, 
certainly,  of  within  10  km. 

Calibration  of  the  3E  Missouri  Region 
Another  approach  to  the  same  question  and  a  benefit  derived  from  the 
BILLIKIW  experiment  is  the  "calibration*’  of  the  southeast  Missotiri 
seismic  region.  Signals  from  the  RTT.T.TKCTJ  experiment  make  it  possible  to 
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Figure  6,  Tinse  differ«ice  hyperbolas,  based  on  calibra¬ 
tion  of  the  southeast  Kissouri  area  by  travel  times 
from  the  BILLIKfc;.  event. 
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determine  the  time  difference  from  a  fixed  point  (the  shot  point)  to  any 
two  pairs  of  the  permanent  seismographic  stations  of  the  region.  Inde¬ 
pendent  of  structsire,  therefore,  with  an  assiimption  only  as  concerns  the 
velocity,  curves  of  equal  time  difference  may  be  plotted.  Figure  6, 
for  instance,  shows  the  time  difference  hyperbolas  for  Holla  and  St. 

Louis  and  for  St.  Louis  and  Bloomington.  The  intersection  of  lines  of 
equal  time  difference  for  two  or  more  pairs  of  stations  permit  the  location 
of  an  epicenter. 

The  time  difference  curves  have  been  used  to  locate  several  earth- 
quoJces  occurring  in  the  southeast  Missouri  seismic  area.  Table  3  compares 
the  time  difference  location,  based  on  the  RTT.T.TKFW  calibration  curves, 

Table  5»  Comparison  of  Least-Square  and  Time-Difference 
Epicenter  Determinations 


Epicenter  Determination 
Least-Square  Time-Difference 


Earthquake 

2  February,  1962 

3  inarch,  1963 

6  Apr  (08^),  1963 


36.493®N,  89.592*W 
36.649*N,  90.026®W 
36.42”  N,  89.72”  W 


36.40”N,  89.63 ”W 
36.65”N,  90.10”W 
36.47”N,  90.00”W 


Sst.  Acc. 
jO.l” 

-tO,2” 


and  the  c<wnputer  determined  epicenters  based  on  knowledge  of  the 
crustal  structure  gained  from  nizmerous  previous  studies  in  tke  region. 

The  BILLIKI2I  based  epicenters  differ  little  from  those  determined  by 
normal  procedures. 

The  use  of  time  difference  curves  has  been  applied  in  the  areas  of 
other  networks.  The  only  distinction  in  the  present  case  is  that  the 
time  differences  are  calibrated  with  reference  to  a  known  source  location. 
In  SE  Missouri  the  knowledge  that  is  at  present  available  of  the  crustal 
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structure  of  the  region  is  such  that  the  time  difference  method  offers 
no  real  adrantage  over  the  usual  P-0  least  square  solution.  'Oie  results 
jnresented  in  table  3  indicate  that  a  calibration  shot  in  a  region  of 
poorly  known  crustal  structure  might  be  used  to  improve  the  accuracy  of 
location  of  a  previously  recorded  event.  In  SE  Missouri  the  time  dif¬ 
ference  curves  also  afford  a  means  for  obtaining  a  rapid  preliminary 
epicenter. 


4.  Other  Studies 

A  comparison  has  also  been  made  of  the  BILLIKEN  data  with  the  re- 
8\ilts  of  regional  crustal  studies.  Ihe  data  themselves  are  not  too  suit¬ 
able  for  use  in  crustal  studies,  for  nowhere  are  there  truly  continuoiis 
profiles  nor  are  the  profiles  reversed.  On  the  other  hand,  the  data  are 
better  than  those  provided  by  any  previous  earthquake  study.  Consequently, 
certadn  comparisons  are  of  advantage. 

Travel  Times 

Figure  7  presents  all  of  the  P  wave  arrivals  on  a  reduced  travel-time 
curve.  Die  times  are  reduced  by  A/6.0  km  per  second.  The  large  circles 
represent  observations  at  sin^e  stations.  A  pair  of  smaller  circles 
joined  by  a  line  segment  represents  an  observation  at  a  spread  of 
seisDKxneters  at  one  of  the  USGS  units.  The  times  of  arrival  of  secondary 
phases  or  of  successive  crests  or  trou^is  after  the  first  breaks  are 
also  plotted. 

The  travel-time  curve  shown  i'.  the  figure  is  based  on  data  from  a 
refraction  survey  and  has  been  proposed  by  Stewart  and  Stauder  (1964) 
as  a  possible  interpretation.  It  assumes  velocity  gradients  in  the 
crust.  As  is  evident,  there  is  general  agreement  between  the  BILLHOI 
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data  and  that  cf  the  refraction  survey.  There  are  some  deviations  about 
the  line  of  first  arrivals.  For  the  most  part,  however,  these  aire  small. 
Three  stations  form  an  exception;  DeSoto,  Missouri,  Parsons,  Tennessee, 
amd  Corinth,  Mississippi.  At  these  stations  the  residuals  araovmt  to  as 
much  as  2  seconds,  Ihese  arrival  times  have  been  carefully  checked.  It 
does  not  seem  possible  to  "explain  away"  these  arrivals  unless  there  is 
scsne  gross  error  in  time  correction  or  in  station  location. 

In  figure  8  the  residuals  of  all  first  arrivals  of  the  P  wave  are 

plotted  as  a  function  of  azimuth.  The  residuals  are  the  observed  arrival 

times  minus  the  theoretical  travel  times  from  the  actual  BILLIKEN  shot 

point  for  the  same  model  used  in  the  previous  portions  of  this  paper, 

that  is,  a  three -layered,  4o  km  crust  with  velocities  of  6,1,  6.4,  and  6.7 

per  sec,  and  with  a  velocity  of  8,15  km/sec.  The  average  of  the 

residuals  is  a  few  tenths  of  a  second  positive.  This  probably  implies  a 

slight  correction  to  the  thickness  of  the  model  or  to  the  P  velocity, 

n 

From  90®  west  to  120®  east  there  is  no  notable  effect  of  azimuth.  To  the 
southeast,  however,  the  two  stations.  Parsons , Tennessee ,  and  Corinth, 
Mississippi,  have  residiials  of  minus  2  seconds,  whereas  to  the  southwest 
at  Searcy  and  Huntsville,  Arkansas,  the  residuals  are  positive  by  almost 
2  seconds.  The  stations  in  these  two  quadrants  au:e  too  few  to  draw  any 
further  conclusions.  However,  the  observations  noted  merit  to  be  borne 
in  mind  in  future  considerations  of  travel  times  in  these  directions. 

Focal  Depth 

The  occurrence  of  the  RTTJ.TKCTJ  shot,  a  Known  surface  event,  also 
presents  an  opportunity  to  investigate  evidence  of  focal  depth  for  earth¬ 
quakes  occxxrring  in  the  same  region.  Figure  9  presents  a  plot  of  the 
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(P  -P  )  tinje  interval  vn.  (S^-P  ).  Th’s  ic,  of  course,  an  application 
of  the  so-called  Thirlaway  nethod.  As  a  first  approach,  a  sinjple  one¬ 
layered  crust  of  4o  km  thickness  and  a  P  velocity  of  6.4  km/sec  was 

S 

considered.  The  lower  of  the  tv/o  diagonal  lines  represents  the  relation 

of  (P  -P  )  vs.  (S  -P  )  for  a  surface  focus  for  this  model.  The  upper 

line  represents  this  same  relation  for  a  focal  depth  of  4o  km.  At 

seven  stations  the  records  of  the  RTT.T.TKF?;  event  permitted  the  reading  of 

P  ,  P  and  3  .  The  corresponding  (P  -P  )  and  (S  -P  )  times  are  represented 
^  B  B  B  ^  S  ^ 

by  the  black  dots  on  the  figure.  As  may  be  seen,  these  points  correspond 
satisfactorily  to  the  surface  focus  line.  For  the  earthquake  of  March  3» 

1963*  P  ,  P  ,  and  S  were  readable  at  another  seven  stations,  and  permit 

’  n*  g*  g  » 

the  plotting  of  the  points  shown  by  the  open  circles.  Data  published  by 

Stauder  and  Bollinger  (1963)  for  the  earthquake  of  the  2rd  of  February, 

1962,  permit  the  plotting  of  the  points  indicated  by  the  crosses.  At  one 

of  these  stations  (Bloomington)  the  point  fell  below  the  surface  focus 

line.  Se-exeimination  of  the  record  showed  that  1.3  seconds  after  the  event 

read  as  the  amplitude  suddenly  increases  fourfold  and  the  record  is 

whited  out  from  that  point  on.  Taking  this  time  as  P  the  point  may  be 

B 

re-plotted  at  the  position  shown.  As  is  usual  in  this  method,  there  is 

a  fair  degree  of  scatter  of  the  data  points  and  the  onsets  of  P  and  S 

B  S 

sure  not  entirely  unambiguous.  However,  the  indication  is  that  the 
earthquake  of  March  3rd  has  a  focus  deep  in  the  crust,  whereas  that  of 
February  2nd  occurred  in  the  upper  part  of  the  crust. 
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